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Automated synthesis from behavioural specifications is an attractive and powerful way of construct- 
ing concurrent systems. Here we focus on the problem of synthesising a membrane system from a 
behavioural specification given in the form of a transition system which specifies the desired state 
space of the system to be constructed. We demonstrate how a Petri net solution to this problem, based 
on the notion of region of a transition system, yields a method of automated synthesis of membrane 
systems from state spaces. 



1 Introduction 

Membrane systems ( lfl9l l20l I2T1 1221 ) are a computational model inspired by the functioning of living 
cells and their architecture and in particular, the way chemical reactions take place in cells divided by 
membranes into compartments. The reactions are abstracted to rules that specify which and how many 
molecules can be produced from given molecules of a certain kind and quantity. As a result, membrane 
systems are essentially multiset rewriting systems. The dynamic aspects of the membrane system model 
including potential behaviour (computations), derive from such evolution rules. 

Petri nets (see, e.g., [5j[23]|25j) are a well-established general model for distributed computation with 
an extensive range of tools and methods for construction, analysis, and verification of concurrent systems. 
Their diverse applications areas include computational and operational foundations for problems and 
issues arising in biology; see for example, lil5l . for a recent comprehensive overview of applications of 
Petri nets in systems biology. 

There are intrinsic similarities between Petri nets and membrane systems. In particular, there exists 
a canonical way of translating membrane systems into Petri nets. This translation is faithful in the 
sense that it relates computation steps at the lowest level and induces in a natural way (sometimes new) 
extensions and interpretations of Petri net structure and behaviour (e.g., inhibitor arcs, localities, and 
maximal concurrency). More details on the relationship between Petri nets and membrane systems can 
be found in, e.g., Il9l[l4l. 



Gabriel Ciobanu (Ed.): Membrane Computing and 
Biologically Inspired Process Calculi (MeCBIC 2012) 
EPTCS 100, 2012, pp. 1-03] doi: 10.4204/E PTCS.100.1I 



© J. Kleijn, M.Koutny, M. Pietkiewicz-Koutny & G.Rozenberg 
This work is licensed under the 
Creative Commons Attribution License. 



2 



Membrane Systems and Petri Net Synthesis 



The strong semantical link between the two models invites to extend where necessary and possible 
existing Petri net techniques and bring them to the domain of membrane systems. An example is the 
process semantics of Petri nets that can help to understand the dynamics and causality in the biological 
evolutions represented by membrane systems (8l[T2l. In this paper, we focus on the synthesis problem, 
that is, the problem of automated construction of a system from a specification of its (observed or desired) 
behaviour. 

Automated synthesis from behavioural specifications is an attractive and powerful way of construct- 
ing correct concurrent systems |[I]|2l|4l|6l|7J[T8l|24l. Here we will re-visit the problem of synthesising a 
Petri net from a behavioural specification given in the form of a transition system. The latter specifies the 
desired state space of the Petri net to be constructed. We will recall a solution to this problem based on 
the notion of region of a transition system. We will then demonstrate how this solution leads to a method 
of automated synthesis of basic membrane systems from state spaces. We also discuss how the proposed 
method could be extended to cope with more complicated kinds of membrane systems. 

2 Preliminaries 

Multisets. A multiset over a finite set X is a function : X — >• N = {0, 1 , 2, . . .}. 6 may be represented 
by listing its elements with repetitions, e.g., 6 = {y,y,z} is such that 9(y) = 2, 9(z) = 1, and 6(x) = 
otherwise. 6 is said to be empty (and denoted by 0) if there are no x such that x G 6 by which we mean 
thatx eX and 9(x) > 1. 

For two multisets 6 and 6' overX, the sum 6 + 6' is the multiset given by (6 + d')(x) = 6(x) + d'(x) 
for all and for k G N the multiset k • 6 is given by (k-6)(x) = k-6 (x) for all x G X . The difference 

6-6' is given by (6 - 6'){x) = max{6{x) - 6'(x),0} for all x G X. We denote 6 < 6' whenever 
6(x) < 6'(x) for all x^X, and 6 < 6' whenever 6 < 6' and 6 ^ 6'. The restriction 6\ z of 6 to a subset 
Z C X is given by 6\z(x) = 6{x) for x G Z, and 6\z{x) = otherwise. The size |0| of 6 is given by 
Exex Q(x)- If / : X — s- Y is a function then f(6) is the multiset over Y such that f(6)(y) = L Ae /-i( v ) 6(x), 
for every jSF. 

Step transition systems. A step transition system over a finite set (of actions) A is a triple TS = 
(Q,£/,qo), where: Q is a set of nodes called states; srf is the set of arcs, each arc being a triple (q, a,q') 
such that q,q' G Q are states and a is a multiset over A; and qo G Q is the initial state. We may write 
q — )■ q' whenever (q, a, q') is an arc, and denote by 



the set of nonempty steps enabled at a state q in TS. We additionally assume that: 

• if q — > q' and q A q" then q' = q" (i.e., TS is deterministic); 

• for every state q G Q, there is a path from go leading to q; 

• for every action a G A, there is an arc g A in TS such that a £ a; and 

• for every state q £ Q,we have g — > iff q = q'. 

Let = {Q,£^,qo) be a step transition system over a set of actions A, and JS' = (Q',£/',q' Q ) be 
a step transition system over a set of actions A'. IS and TS' are isomorphic if there are two bijections, 
: A — > A' and V : Q — ^ Q', such that v(qo) = q' and, for all states q,q' G Q and multisets a over A: 



tsSteps q = {a | a ^ A 3q f : q ^ q'} 
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We denote this by TS ~ , v TS' or TS ~ PS'. 



Petri nets. A Place/Transition net (or PT-net) is specified as a tuple PP = (P, T, W,Mq), where: P and 
T are finite disjoint sets of respectively places and transitions; W : (P x P) U (P x T) — >• N is the arc 
weight function; and Mo : P — > N is the initial marking (in general, any multiset of places is a marking). 
We assume that, for each transition t, there is at least one place p such that W(p,t) > 0. In diagrams, 
such as that in Figured] places are drawn as circles, and transitions as boxes. If W(x,y) > 1, then (x,y) is 
an arc leading from x to y. An arc is annotated with its weight if the latter is greater than one. A marking 
M is represented by drawing in each place p exactly M{p) tokens (small black dots). 




Figure 1 : A PT-net. 

A step U of PT is a multiset of transitions. Its pre-multiset and post-multiset of places, *U and U', 
are respectively given by 

'U{p) = 1 £U(t)-W(p,t) and U'(p) = 1 £U(t)-W(t,p), 

teU tell 

for each place p. For the PT-net in Figure[T]we have: 

'{t?\t?\t T »} = {p?,p?,pf} and K",t?\tf }' = {pi,pi,pi,pf,pi,p^}. 

We distinguish two basic modes of execution of PT-nets. To start with, a step of transitions U is 
free-enabled at a marking M if 'U < M. We denote this by M[U)f ree , and then say that a free-enabled U 
is max-enabled at M if U cannot be extended by a transition to yield a step which is free-enabled at M, 
i.e., there is no t G T such that M[U + {t})f ree . We denote this by M[U) max . In other words, U is free- 
enabled at M if in each place there are sufficiently many tokens for the specified multiple occurrence of 
each of its transitions. Maximal concurrency (max-enabledness) means that extending U would demand 
more tokens than M supplies. For the PT-net in Figure Q] we have that, at the given marking Mo, the step 
{t j 12 , tj 21 } is free-enabled but not max-enabled, and {t j 1 1 , t ^ 12 , tj 21 , tj 22 } is max-enabled. 

For each mode of execution m G {free, max}, a step U which is m-enabled at a marking M can be 
m-executed leading to the marking M' given by M' = M — *U + U*. We denote this by M[U) m M'. For 
the PT-net in Figure Q] we have 

M [{t j 12 , tl 21 }}free {Pi , Pi , P2 > P2 . P2 > P2> P3 } ■ 
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Petri nets with localities. PT-nets are a general model of concurrent computation. To capture the com- 
partmentisation of membrane systems, lfl2l adds explicit localities to transitions. Though not necessary 
from a modelling point of view, we associate in this paper — only for notational convenience — also 
each place with a locality. 

A Yl-net with localities (or PTL-net) is a tuple PTL = (P, T,W,£,M ) such that (P, T,W,M ) is a PT- 
net, and £ is a location mapping for the transitions and places. Whenever £(x) = £{£), we call x and z 
co-located. In diagrams, nodes representing co-located transitions and/or places will be shaded in the 
same way, as shown in Figure El 




Figure 2: A PTL-net corresponding to a basic membrane system, where £{x\) = i, for each node of the 
form x\. Note that, e.g., transitions tj" , t j 12 and t^ 13 are co-located. 

Co-locating transitions leads to one more way of enabling for steps of transitions. We say that a 
step U of PTL is Imax-enabled at a marking M if M[U)f ree and U cannot be extended by a transition 
co-located with a transition in U to yield a step which is free-enabled at M; i.e., there is no t € T such 
that £(t) G £(U) and M[U + {t})f ree . We denote this by M[U)i max , and then denote the lmax-execution of 
U by M\U)i max M' , where M' = M — *U + U° . Note that locally maximal (Imax) concurrency is similar to 
maximal concurrency, but now only active localitie^] cannot execute further transitions. For the PTL-net 
in Figure|2]we have that {tj 11 ,tj 12 } is lmax-enabled at the given marking, but {t^ 11 } is not. 

Let m G {free, max, Imax} be a mode of execution of a PTL-net PTL. Then an m-step sequence is 
a finite sequence of m-executions starting from the initial marking, and an m-reachable marking is any 
marking resulting from the execution of such a sequence. Moreover, the m-concurrent reachability graph 
of PTL is the step transition system: 

CRG m (PTL) = ( [M ) m , {(M,U,M') | M £ [M ) m A M[U) m M'} , M ) , 

where [M ) m is the set of all m-reachable markings which are the nodes of the graph; M is the initial 
node; and the arcs between the nodes are labelled by m-executed steps of transitions. Concurrent reacha- 
bility graphs provide complete representations of the dynamic behaviour of PTL-nets evolving according 
to the chosen mode of execution. 



By active localities of a step U we mean the localities of transitions present in U. 
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Figure 3: A membrane structure (m = 3) and its compartments with 1 being the root node, (1,2) G }l 
and 1 = parent(3). 

Membrane structures. A membrane structure /J. (of degree m > 1) is given by a rooted tree with m 
nodes identified with the integers 1, .. . ,m. We will write G jj, or i = parent(j) to indicate that there 
is an edge from i (parent) to j (child) in the tree of pL, and i G /I means that i is a node of pt. The nodes of 
a membrane structure represent nested membranes which in turn determine compartments (compartment 
i is enclosed by membrane i and lies in-between i and its children, if any), as shown in Figure [3] 

We will say that a PTL-net PTL = (P,T,W,£,Mq) is spanned over the membrane structure \l if £ : 
P U T -> ;U and the following hold, for all /? G P and t G T: 

• if W(p,f) > then £{p) = £(t); and 

• if W(t,p) > then £{p) = £{t) or (£(p),£(t)) G M or (£(t),£(p)) G At. 

The PTL-net of Figure |2] is spanned over the membrane structure depicted in Figure [3] 

Basic membrane systems. Let V be a finite alphabet of names of objects (or molecules) and let fx be 
a membrane structure of degree m. A basic membrane system (over V and pt) is a tuple 

SMS=(V,At,w?,...,w^/? 1 ,...,/? m ) 

such that, for every membrane i, is a multiset of objects from V, and R{ is a finite set of evolution 
rules associated with membrane (compartment) i. Each evolution rule r G /?,• is of the form r : lhs r — >• r/is r , 
where //z/ (the left hand side of r) is a nonempty multiset over V, and r/2/ (the right hand side of r) is a 
multiset over 

V U {a out \ aeV}U {ai nj \ a G V and G jJ.} . 

Here a symbol a;„ represents an object a that is sent to a child node (compartment) j and a out means that 
a is sent to the parent node. If i is the root of p, then no indexed object of the form a out belongs to rhs r . 
A configuration of BMS is a tuple 

C= (wi,...,w m ) 

of multisets of objects, and Co = (wj, . . . , w° ) is the initial configuration. Figure @] shows a basic mem- 
brane system over the membrane structure depicted in Figure [3] 

A membrane system evolves from configuration to configuration as a consequence of the application 
of evolution rules. There are different execution modes ranging from fully synchronous — as many 
applications of rules as possible — to sequential — a single application of a rule at a time. Here, similarly 
as in the case of PTL-nets, we distinguish three modes, all based on the notion of a vector multi-rule. 
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{b} -4 {a} 

{a} -4 {b,Ci„ 2 ,a in3 } 

{b} -> {c,a ( „ 3 } 



1 




2 




{a,Z?,c,c} 



r2i : {a,c} — > {£>} 
r 22 : {£} -)• {a} 



Figure 4: Basic membrane system BMSq. 



A vector multi-rule of BM5 is a tuple r = (ri,...,r m ) where, for each membrane i of il, r, is a 
multiset of rules from /?,. For such a vector multi-rule, we denote by lhs\ the multiset 



in which all objects in the left hand sides of the rules in r, are accumulated, and by rhs] the multiset 



of all (indexed) objects in the right hand sides. The first multiset specifies how many objects are needed 
in each compartment for the simultaneous execution of all the instances of evolution rules in r. 
A vector multi-rule r of BMS is 

• free-enabled at a configuration C if Ihs] < w,-, for each i. 
Moreover, a free-enabled vector multi-rule r = (ri,... ,r m ) is: 

• max-enabled if no r, can be extended to a vector multi-rule which is free-enabled at C; and 

• Imax-enabled if no nonempty r, can be extended to a vector multi-rule which is free-enabled at C. 
For example, in Figure |4j 

• (0,0, {r3i }) is not free-enabled; 

• ({rn , r^}, 0,0) is lmax-enabled but not max-enabled; and 

• ({ri 1 ,ri 2 },{r2i,r 2 2},0) is max-enabled. 

If r is free-enabled (free) at a configuration C, then C has in each membrane i enough copies of 
objects for the application of the multiset of evolution rules r,-. Maximal concurrency (max) requires that 
adding any extra rule makes r demand more objects than C can provide. Locally maximal concurrency 
(Imax) is similar but in this case only those compartments which have rules in r cannot enable any more 
rules; in other words, each compartment either uses no rule, or uses a maximal multiset of rules. 



£ rj(r)-rhs r 

reRi 
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The effect of the rules is independent of the mode of execution m € {free, max, Imax}. A vector 
multi-rule r which is m-enabled at C can m-evolve to a configuration C = (w\ w' m ) such that, for each 
i and object a: 

w'i(a) = Wj (a) - Ihs] (a) + rhs] (a) + rhs x parent[i) {a ini ) + £ rhs* (a out ) 

i=parent(j) 

where rhs parent ^ = if i is the root of jX. We denote this by C — ^ m C . Moreover, an m-computation is 
a finite sequence of m-evolutions starting from the initial configuration; any configuration which can be 
obtained through such a computation is called m-reachable. For the basic membrane system depicted in 
Figure @] we have, for example: 

({rii,r 12 },0,0) ff , , , , , , {0,{r 2 i,r 2 2},0) , f r , -, f n 

Let m S {free, max, Imax} be a mode of execution of a basic membrane system BMS. Then the 
m-concurrent reachability graph of BMS is given by: 

CRG m (BMS) = ( [C ) m , {(C,ri + . . . + r m ,C') | C G [C ) m A C i^^H m C} , C ) , 

where [Co) m is the set of all m-reachable configurations which are the nodes of the graph; Co is the initial 
node; and the arcs between the nodes are labelled by multisets of evolution rules involved in the m- 
executed vector multi-rulesH Similarly as in the case of PTL-nets, concurrent reachability graphs capture 
completely the dynamic behaviour of basic membrane system evolving according to the chosen mode of 
execution. 



3 Membrane Systems and Petri Nets 



There is a natural way of translating a basic membrane system BMS = (V, jl , w\ , . . . , w% , Ri , . . . , R m ) over 
a membrane structure /I into a behaviourally equivalent PTL-net PTL(BMS) = {P,T,W,1,Mq) spanned 
over the same membrane structure. In the constructed net, places represent objects present inside com- 
partments, and transitions represent evolution rules. Both places and transitions are associated with 
membranes and this information is represented by the location mapping. 

The constructed PTL-net PTL(BMS) has a separate place p" with £(p") = j, for each object a and 
membrane j, and a separate transition tj" with l{t*) = i, for each rule r in compartment i. 

The initial marking inserts w Aa) tokens into each place p°. The connectivity between transition 
t = t- and place p = p" is given by: 



as well as: 



W(p,t) 



W(t,p) 



Ihs' (a) if i = j 
otherwise 



rhs r (a) 
rhs r {a out ) 
rhs r (a 




if i = 

if J 

if i-- 



J 



parent (i) 
, . parent (j) 
otherwise . 



2 Though it may be that rules from different membranes are the same in terms of the multisets defining their left hand and 
right hand sides, we assume here that evolution rules associated with different membranes can be distinguished, e.g., by giving 
them each their own name (an injective label). 
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Figure |2] shows the result of the above translation for the basic membrane system in Figure |4] Note that 
it immediately follows from the construction that the PTL-net PTL(BMS) is spanned over jj.. 

The PTL-net PTL(BMS) provides a. faithful representation of the behaviour of the basic membrane 
system BMS. To capture this very close relationship, we define two bijective mappings, v and p, which 
allow us to move between BMS and PTL(BMS): 

• for every marking M of PTL(BMS), v(M) = (wi,. . ■ ,w m ) is the configuration of BMS given by 
Wi(a) = M (pf ), for every object a and every i. 

• for every step U of PTL(BMS), p(U) = (r\,...,r m ) is the vector multi-rule of BMS given by 
r,(r) = U (tj ), for every rule r G /?, and every i. 

It is then possible to establish a direct relationship between (the operation of) the original membrane 
system and the PTL-net resulting from the above translation at the system level: 

C-^ m C v-\C)[p-\r)) m v-\C) 

M[U) m M> V(M) ^ m v(M') 

for all modes of execution m G {free, max, Imax}, configurations C of BMS and markings M of PTL{BMS). 
Together with v(Mq) = Co, this result means that the m-step sequences of PTL{BMS) faithfully represent 
m-computations of BMS, and the same applies to markings and configurations. Crucially, we obtain 

Theorem 1 For each m G {free, max, Imax}, 

CRG m (PTL(BMS) ) CRG m (BMS) , 

where the mapping v is defined as above, and 0(tf) = r,/or every transition t\ of PTL{BMS). 

The above theorem captures the very tight behavioural correspondence between BMS and PTL{BMS), 
allowing to apply analytical techniques developed for Petri nets in the analysis of membrane systems. 
For example, one can employ the invariant analysis based on linear algebra (26), or use the causality se- 
mantics approach of Petri nets based on occurrence nets, as first outlined in [12]|. In this paper, we show 
how techniques used to synthesise Petri nets could be employed in order to construct basic membrane 
systems from their intended behaviours as represented by step transition systems. First, however, we 
provide a translation from PTL-nets spanned over membrane structures to basic membrane systems. 

Let PTL = (P,T,W,£,Mo) be a PTL-net spanned over a membrane structure pL. For such a PTL-net, 
we construct the corresponding basic membrane system BMS(PTL) over ji. in the following way: 

• P is the set of objects; 

• the initial configuration is v'(Mo) where, for every marking M of PTL, 

V '(M) = (M\ e - l{l)nP , . . . ,M\t-\m)np) ! 

• each transition t G T with t* = {p 1 ,. . . ,p } has a corresponding evolution rule <j)'(t) of the form 
t : *t — > {a\ , . . . ,a,k} where, for i = l,...,k, 







if £(/?'') 


= l{t) 


<k = < 


Pout 


if lip 1 ) 


= parent (£(t)) 






if 4(f) 


= parent(£(p 1 )) 
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• for each membrane i G pt, the set of evolution rules is given by R, = {(j>'(t) \ £{t) = i}. 
Again, the translation results in a very close behavioural correspondence. 

Theorem 2 For each m G {free, max, Imax}, 



CRG m (PTL) ~ ¥y CRG m (BMS(PTL) ) 



where the mappings (f)' and v' are defined as above. 

It follows from Theorems Q] and [2] that the problem of synthesis of basic membrane systems from 
step transition systems is equivalent to the problem of synthesis of PTL-nets spanned over membrane 
structures. It therefore suffices to solve the latter, and in the next section we describe a solution based on 
the notion of a region of a step transition system. 

4 Synthesising nets corresponding to membrane systems 

The Petri net synthesis problem we consider is formulated as follows. 

Problem 1 Given are a finite set T, a membrane structure pi, a mapping £ : T — >■ pt, m G {free, max, Imax}, 
and TS = {Q,si ,qo) which is a finite step transition system over T. 

Construct a PJh-net PTL = (P, T,£,M ) spanned over jj. such that CRG m (PTL) ~ TS, and £ is an exten- 
sion of the mapping defined for T. 

As demonstrated in [4], synthesis problems like Problem Q] can be solved using techniques coming 
from the theory of regions of transition systems (see, e.g., HI 13 [HI)- Intuitively, a region represents a 
single place in a hypothetical net generating the given transition system. Regions are used both to check 
whether a net satisfying the conditions can be constructed and, if the answer turns out to be positive, to 
construct such net. 

In this particular case, a region of the step transition system TS consists of three mappings 



Here co(a) = Y*teT a if) ■ co(t) and similarly i(a) = Y*teT a { t ) " m a re gi° n of the form (0 repre- 
senting a place p, o{q) is the number of tokens in p in the marking corresponding to the node q, co(t) 
represents the weight of the arc from p to transition t, and i(t) represents the weight of the arc from t 
to p. It is then natural to require in (O that p contains enough tokens not to block a step a executed at q, 
and also to ensure that the number of tokens in p before and after executing a is consistent with the total 
arc weight of the step a in relation to p. 

In the case of Problem [T] one also needs to take into account the fact that the target PTL-net must 
be spanned over pt. This imposes additional constraints on allowed regions (places) and the location 
mapping I. We call a region reg as in (O with a location l{reg) G \l compatible with the membrane 
structure \i if the following hold, for every t S T: 

• if co(t) > then £(t) = £(reg); and 




(2) 



such that, for every arc q — > q' °f TS, 



o(q)>G)(a) and o{q) = c(q) - w(a) + i(a) . 



(3) 



• if i(t) > then £(t) = £{reg) or (£(t),£{reg)) G /I or (£(reg),£{t)) G pt. 
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The set of all such regions will be denoted by P^. Note that if reg is such that co(t) > 0, for at least one 
(GT, then l(reg) is uniquely determined; otherwise we always choose l{reg) to be the membrane which 
is higher up in the tree structure of 11 than any other suitable candidate. As a result, we can leave &(reg) 
implicit. 

Finally, Problem Q] should be feasible in the sense that the transition system TS can be realised by a 
suitable net. There are two necessary and sufficient conditions for realisability (see fflfTTI): 

• state separation: for every pair of distinct states of the transition system there is a region (a marked 
place) distinguishing between them; and 

• forward closure: there are sufficiently many places defined by regions of the transition system to 
disallow steps not present in the transition system. 

First we describe how all places can be found that potentially provide a solution to Problem [Q in other 
words, all the regions Q of the transition system TS which are compatible with 11. 

Finding compatible regions. Let T, [i, £:T fi and TS = (Q,£^,qo) be as in Problem [TJ Assume 
that Q = {qo, ...,q/ r } and T = {h, . . . ,?„}. We use three vectors of non-negative variables: 

x = x ...x h y = y 1 ...y„ z = z\...z n - 

We also denote p = xyz and define a homogeneous linear system 

^:( X '- a ' Z . , for all q x A q, in TS 

where a ■ z denotes a{t\) -z\ H V cc(t n ) ■ z. n and similarly for a ■ (y — z). 

The regions (0 of TS are then determined by the integer solutions p of the system & assuming that, 
for < i < h and 1 < j < n, 

o(qi) = Xi i(tj) = yj (0(tj) = Zj 

The set of rational solutions of forms a polyhedral cone in Q /l + 2,I + 1 . As described in (3), one 
can effectively compute finitely many integer generating rays p 1 , . . . ,p of this cone such that any in- 
teger solution p of can be expressed as a linear combination of the rays with non-negative rational 
coefficients: 

k 

i=i 

Such rays p' are fixed and (some of them) turned into net places if Problem Q] has a solution. More 
precisely, if p' is included in the constructed net, then 

M (p')=4 W(p l ,t i )=z' i W(t i ,p')=y l i , (4) 

where Mo is the initial marking of the target net, and ti&T. 

Clearly, not all such rays can be considered for the inclusion in the net being constructed, as the 
corresponding regions have to be compatible with pL. We therefore ensure through a simple check that 
the generating rays p 1 , . . . ,p* are compatible with jx, deleting in the process those which are not. Note 
that any p € is a non-negative linear combination of rays compatible with ii. 

Having found the generating rays compatible with jj., we proceed to check whether Problem Q] has 
any solutions at all. 
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Checking state separation. Let TS = (Q, £?,qo) be as in ProblemQ] We take in turn each pair of dis- 
tinct states, qt and qj, of Q and decide whether there exists p = (a, i, (o) G with coefficients c\ , . . . ,ct 
such that c(qi) = Xi ^ xj = o(qj). Since the latter is equivalent to 



1=1 1=1 

one can simply check whether there exists at least one p' (called a witness (H) such that x\ ^ x!j. 

Checking forward closure. Again, let TS = (Q, qo) be as in ProblemQ] and m G {free, max, Imax}. 
First, we take in turn each state qi of Q, and calculate the set of region enabled steps, denoted by 
regSteps q . . Intuitively, region enabled steps are those that cannot be disabled (or blocked) by compatible 
regions. 

To build regSteps q . one only needs to consider nonempty steps a with | a | < m ■ Max, where Max is 
the maximum size of steps labelling arcs in TS, and m is the number of membranes of }X. The reason is 
that, for each membrane i G jit there exists a compatible region (a, i, <o) G (called a witness) such that 
g{Q) = {Max} and, for every t G T , 



(o(t) =l(t) 



1 if £(t) = i 
otherwise . 



Taken together, all such regions block any step a with | a \ > m ■ Max. 

For each nonempty step a with \a\ < m-Max it is the case that a ^ regSteps q . iff for some p G P^ 
with coefficients ci, . . . ,Ck we have Xi < a • z. Since the latter is equivalent to 



k 

Y^ci ■ (x\ - a -z') < 

1=1 



one simply checks whether there exists at least one p' (again called a witness) such that x\ — a ■ z 1 < 0. 

Having determined the region enabled steps, in order to establish forward closure we need to verify 
that, for every state q G Q, 



tsSteps 



regSteps q if m =free 

{a G regSteps q \ —3t G T : a + {t} G regSteps q } if m = max 

{a G regSteps q \ -Bt G T : a + {?} G regSteps q A G ^(a)} if m = /max . 



Constructing the solution net. If the above checks for the feasibility of Problem Q] are successful, one 
can construct a solution PTL-net spanned over pt by taking all the witness rays and regions, and treating 
them as places in the way indicated in (0]). The resulting net PTL satisfies 

CRG m {PTL) ~ TS . 



5 Concluding remarks 

We have described how one can adapt a solution to the Petri net synthesis problem based on regions of 
step transition systems, so that the resulting method can be used to construct basic membrane systems 



12 



Membrane Systems and Petri Net Synthesis 



with a specific behaviour. Moreover, there are other synthesis results developed for Petri nets which can 
be employed to extend the proposed solution in several directions, two of which are briefly mentioned 
below. 

In Problem Q] it is assumed that the association of transitions with membranes is given. This can 
be relaxed and one can aim at synthesising membrane systems without such an association, or even 
without being given a membrane structure (in such a case, the synthesis procedure should construct a 
membrane structure as well). For such a modification, there already exist results which can be used to 
develop a solution. More precisely, the method of 'discovering' localities in iTTTl works for PTL-nets 
with localised conflicts (where transitions which share an input place are co-located). Since all PTL-nets 
spanned over membrane structures have localised conflicts, the result in [ 17 ] can be adapted to work for 
basic membrane systems. 

Evolution rules of membrane systems are often equipped with promoters and inhibitors. Both fea- 
tures have direct counterparts in Petri nets in the form of activator and inhibitor arcs, and suitable transla- 
tions between membrane systems and Petri nets can be developed as described in (UdH. Moreover, the 
synthesis technique based on regions of step transition systems works also for PTL-nets extended with 
activator and inhibitor arcs iPToTl . In fact, there is a general setting of so-called T-nets and corresponding 
T-regions HI SI. Here the parameter % is a general and convenient way of capturing different types of 
connections (arcs and their combinations) between places and transitions, removing the need to re-state 
and re-prove the key results every time a new kind arcs is introduced. Note that the recently introduced 
SET-nets lfT3l [T4l (with qualitative rather than quantitative resource management) and set membrane 
systems iflOl can be treated within the general theory of T-net synthesis based on regions of transition 
systems fTTI . 
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